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The effective evolution of an inhomogeneous universe model in Einstein's theory of gravitation 
may be described in terms of spatially averaged scalar variables. This evolution can be modeled by 
solutions of a set of Friedmann equations for an effective scale factor, with matter and backreaction 
source terms, where the latter can be represented by a minimally coupled scalar field ('morphon 
field') . We review the basic steps of a description of backreaction effects in relativistic cosmology that 
lead to refurnishing the standard cosmological equations, but also lay down a number of unresolved 
issues in connection with their interpretation within observational cosmology. 
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I. GENERAL THOUGHTS: THE STANDARD 
MODEL, THE AVERAGING PROBLEM AND 
KEY INSIGHTS 

A. The standard model of cosmology 

The standard model of cosmology does not, like the 
standard model of particle physics, enjoy appreciable 
generality; it is based on the simplest conceivable class 
of (homogeneous-isotropic) solutions of Einstein's laws of 
gravitation. It is clear that the inhomogeneous properties 
of the Universe cannot be described by such a strong ide- 
alization. The key issue is whether it can be described so 
on average, and this is the subject of considerable debate 
and controversy in the recent literature. If the standard 
model indeed describes the averaged model, we have to 
conjecture that backreaction effects, being the main sub- 
ject of this contribution, are negligible. We are striving 
to discuss most of the related aspects of this debate. 



1. Dark Energy and Dark Matter 

In the standard model of cosmology one has to con- 
jecture the existence of two constituents, if observational 
constraints are met, that both have yet unknown ori- 
gin: first, a dominant repulsive component is thought 
to exist that can be modeled either by a positive cos- 
mological constant or a scalar field, e.g. a so-called 
quintessence field, possibly deriving either from higher- 
order Ricci curvature Lagrangians [3j|,_|29| or string- 
motivated low-energy effective actions In any case 
we should not expect a fundamental scalar field to ex- 
ist in nature, at least not one that can be viewed as a 
natural candidate for the relevant effects needed. This 
expectation is also hampered by the well-known viola- 
tion of energy conditions of a quintessence field that is 
able to produce late-time volume acceleration of the Uni- 
verse. Rather, a scalar field would likely be an effective 
one, either stemming from higher-order gravity terms, or 
effective terms as remnants from higher dimensions that 
are compactificd or even non-compactified as in brane 
world cosmologies [6(| ■ As we shall learn below, already 
classical general relativity allows to identify effective ge- 
ometrical terms, simply resulting from inhomogeneities, 
with an effective scalar field component, the morphon 
field [26[ , a good example of William of Ockham's razor. 
If we restrict our attention to cosmology and the fitting 
of extra terms from various different theories to observa- 
tional data, then those extra terms may also be mapped 
into morphon fields with different but unambiguously de- 
fined physical consequences. A review of the status and 
properties of such models can be found in [3||, see also 

MM- 

Besides this Dark Energy, there is, secondly, a non- 
baryonic component that should considerably exceed the 
contribution by luminous and dark baryons and massive 
neutrinos. This Dark Matter is thought to be provided by 



exotic forms of matter, not yet detected in experiments. 
According to the concordance model (57| , the former com- 
prises almost 3/4 and the latter almost 1/4 of the total 
source of Friedmann's equations, up to a few percent that 
have to be attributed to baryonic matter and neutrinos 
(in the matter-dominated era). The intriguing question 
is whether an explanation of these missing components 
is the task of particle physicists, or whether the cosmo- 
logical model is built on oversimplified priors. 

2. Uncharted territory beyond the standard model 

The concordance model is encircled by a large set 
of observational data that are, however, orthogonal 
only within the predefined solution space of a FRW 
(Friedmann-Robertson-Walker) cosmology. This solu- 
tion space is two-dimensional, for Friedmann's expan- 
sion law derives from Hamilton's constraint restricted to 
locally isotropic and hence (by Schur's Lemma) homoge- 
neous distributions of matter and curvature, 

n m + n k + n A = 1 , (i) 

where the standard cosmological parameters are global 
and iconized by the cosmic triangle [7j, 

_ SttGqh . -k _ _A_ 

m 3ff 2 ' fc a 2 ff 2 ' A 3i? 2 ' 1 j 

Quit) is the homogeneous matter density, H(t) := a/ a 
Hubble's function with the scale factor a(t), k a posi- 
tive, negative or vanishing constant related to the three 
elementary constant-curvature geometries, and A is the 
cosmological constant. 

We shall learn below that an extended solution space 
of an averaged inhomogeneous universe model is already 
three-dimensional, when we include inhomogeneities of 
matter and geometry. Hence, such more realistic models 
enjoy more parameter freedom, but at the same time de- 
fine these parameters physically in terms of volume aver- 
ages. How can we be sure that fitting an idealized model, 
that ignores inhomogeneities, to observational data is not 
'epicyclic', especially if the model enters as a prior into 
the process of interpreting the data? Confronting ob- 
servers with the wider class of averaged cosmologies al- 
lows them to draw their data points within a cube of 
possible solutions and to differentiate the relevant ob- 
servational scales reflected by these data; if we 'force' 
them to draw the data points into the plane of the FRW 
solutions on every scale, then they conclude that there 
are 'dark' components. Thus, we have to exclude that 
they may have missed something in the projection and we 
have to clarify whether the ignorance of scale-dependence 
of observables in the standard model does not mislead 
their interpretation. Both issues are equally important to 
judge the viability of the standard model in observational 
cosmology: the first is the question of how backreaction 
quantitatively affects the standard cosmological parame- 
ters, and the second is the comparison of data taken on 
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small scales (e.g. on cluster scales) and data taken on 
large scales (e.g. CMB; high-redshift supernovae). Both 
additional 'degrees of freedom' in interpreting observa- 
tional data are interlocked in the sense that backreaction 
effects may alter the evolution history of cosmological pa- 
rameters. A comparison of data taken on different spatial 
scales has therefore also to be subjected to a critical as- 
sessment of data that are taken at different times of the 
cosmic history: with backreaction at work, the simple 
time-scaling of parameters in a FRW cosmology is also 
lost. 

Before we enter the physics of backreaction that is easy 
to understand, we have to first probe some more critical 
territory in the following subsection. 



B. Different 'directions' of backreaction 

The notion of averaging in cosmology is tied to space- 
plus-time thinking. Despite the success of general co- 
variance in the four-dimensional formulation of classical 
relativity, the cosmologist's way of conceiving the Uni- 
verse is evolutionary. This breaking of general covari- 
ance is in itself an obstacle to appreciating the proper 
status of cosmological equations. The standard model of 
cosmology is employed with the implicit understanding 
that there is a global spatial frame of reference that, if 
mapped to the highly isotropic Cosmic Microwave Back- 
ground, is elevated to a physical frame rather than a 
particular choice of a mathematical slicing of spacctimc. 
Restricting attention to an irrotational cosmic continuum 
of dust (that we shall mostly retain throughout this con- 
tribution), the best we can say is that all elements of the 
cosmic continuum defined by the homogeneous distribu- 
tion of matter are in free fall within that spacetime, and 
therefore are preferred relative to eventually accelerating 
observers with respect to this frame of reference. Those 
preferred observers are called fundamental. Exploiting 
the diffeomorphism degrees of freedom we can write the 
FRW cosmology in contrived ways, so that nobody would 
realize it as such. This poi nt is raised as a criticism of 
an averaging framework [5fJ , as if this problem were not 
there in the standard model of cosmology. Again, the 
'natural' choice for the matter model 'irrotational dust' 
is a collection of free-falling continuum elements, now for 
an inhomogeneous continuum. 

Let us consider the 4-metric form to be 



g 



-dr 



3 g ; 3 g = gab dx a ® dx b 



(3) 



where latin indices run through 1 • • • 3 and X a are lo- 
cal (Gaussian normal) coordinates. Evolving the first 
fundamental form 3 g of the spatial hypersurfaces along 
d/dt —: dt defines their second fundamental form 

3 K = K ab dX a ® dX b ; K ab := -~d t9ab , (4) 

with the extrinsic curvature components K ab . Such a co- 
moving ( synchronous ) slicing of spacetime may be consid- 



ered 'natural', but it may also be questioned. However, 
to dismiss its physical relevance due to the fact that shell- 
crossing singularities arise is shortsighted. It is a problem 
of the matter model in the first place. A comoving (La- 
grangian) frame helps to access nonlinear stages of struc- 
ture evolution, as is well-exemplified in Newtonian mod- 
els of structure formation, where the 'slicing problem' is 
absent. Those nonlinear stages inevitably include the de- 
velopment of singularities, provided we do not improve on 
the matter model to include forces that counteract grav- 
ity (like velo city dispersion) in order to regularize such 
singularities 23J . If a chosen slicing appears to be better 
suited, because it does not run into singularities, then 
one should rather ask the question whether the evolution 
of variables is restricted to a singular-free regime just be- 
cause inhomogeneities are not allowed to enter nonlinear 
stages of structure evolution. An example for this is per- 
turbation theory formulated e.g. in longitudinal gauge, 
where the variables are 'gauge-fixed' to a (up to a given 
time-dependent scale factor) non-evolving background. 

However, the slicing problem is not off the table. There 
exist strategies to consolidate the notion of an effective 
spatial slicing that would minimize frame fluctuations be- 
ing attributed to the diffeomorphism degrees of freedom 
in an inhomogeneous model. Such, more involved, strate- 
gies relate to the intrinsic direction of backreaction that 
we put into perspective below. 



1. Extrinsic (kinematical) and intrinsic backreaction 

Having chosen a foliation of spacetime implies that we 
can speak of two 'directions': one being extrinsic in the 
direction of the extrinsic curvature K ab of the embedding 
of the hypersurface into spacetime (e.g. parametrized by 
time), the other being intrinsic in the direction of the 
Ricci tensor R ab of the three-dimensional spatial hyper- 
surfaces parametrized by a scaling parameter (let it be 
the geodesic radius of a randomly placed geodesic ball). 
Consequently, we may speak of two 'directions' of back- 
reaction: inhomogeneities in extrinsic curvature and in 
intrinsic curvature. The former is of kinematical nature, 
since we may interpret the extrinsic curvature actively 
through the expansion tensor <d ab := —K ab , and intro- 
duce a split into its kinematical parts: Q a b — \g a b < d+< J ab, 
with the rate of expansion = C C , the shear tensor o~ ab , 
and the rate of shear a 2 := l/2a ab a ab ; note that vorticity 
and acceleration are absent for dust in the present flow- 
orthogonal foliation. The latter addresses the so-called 
fitting problem [38l.l42l]. [68| . i.e. the question whether we 
could find an effective constant-curvature geometry that 
replaces the inhomogeneous hypersurface. An answer to 
this question has to deal with the problem of 'averag- 
ing' the tensorial geometry for which several different 
strategies are conceivable. Some of those strategies do 
not distin guis h between extrinsic and intrinsic averaging 
(e.g. [sH, 34 1, and other references in 39]). One method 



has recently obtained a strong position in the context of 
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Perelman's work (e.g. [70|) on the Ricci-Hamilton flow 
related to the recent proof of Poincare's conjecture, and 
implied progress on Thurston's geometrization program 
to cut a Riemannian manifold into 'nice pieces' of eight 
elementary geometries. This method we briefly scetch 
now. 



Smoothing the geometry 



C. The origin of kinematical backreaction and the 
physics behind it 

Let us now concentrate on the question, why there 
must be backreaction at work, restricting attention to 
kinematical backreaction as defined above. After we have 
understood the reasons behind backreaction effects in 
general terms, i.e. without resorting to restrictions of 
spatial symmetry or approximations of evolution mod- 
els, the very question of their relevance is better defined. 



Employing the Ricci-Hamilton flow [47|, [48j|, an 'av- 
eraging' of geometry can be put into practice by a rescal- 
ing of the metric tensor. A general scaling flow is de- 
scribed by Petersen's equations [7l| that we may imple- 
ment through a 2+1 setting by evolving the boundary 
of a geodesic ball in a three-dimensional cosmological 
hypersurface in radial directions, thus exploring the Rie- 
mannian manifold passively. Upon linearizing the gen- 
eral scaling flow, e.g. in normal geodesic coordinates, we 
obtain a scaling equation for the metric along radial di- 
rections [l9j : up to tangential geometrical terms on the 
boundary we obtain: 



-2K ab (r Q )[r - r ] 



(5) 



i.e. the metric scales in the direction of its Ricci tensor 
much in the same way as it is deformed in the direction 
of the extrinsic curvature by the Einstein flow. If we now 
implement the active (geometrically Lagrangian) point 
of view of deforming the metric by the same flow along a 
Lagrangian vector field d/drg while holding the geodesic 
radius ro fixed, we are able to smooth the metric in a 
controlled way. Depending on our choice of normaliza- 
tion of the flow, we may preserve the mass content inside 
the geodesic ball while smoothing the metric. Such a 
rescaling transforms average properties on given hyper- 
surfaces from their values in the inhomogeneous geometry 
(the actual space section) to their values on a constant- 
curvature geometry (the fitting template for the space 
section): they are renormalized resulting in additional 
backreaction effects due to the difference of the two vol- 
umes (the Riemannian volume of the actual space section 
and the constant-curvature volume) - the volume effect, 
and also curvature backreaction terms that involve av- 
eraged invariants of the Ricci tensor. For details and 
references see and for small overviews [2(| and [2l| . 
In such a setting the role of lapse and shift functions (i.e. 
the slicing problem) can also be controlled by employing 
the recent results of Perelman [3l[ . 

We now come to some crucial points of understanding 
the physics behind backreaction. In order not to think 
of any exotic mechanism, the historical use of the notion 
'model with backreaction' should simply be replaced by 
'more realistic model'. 



1. An incomplete message to particle physicists 

Employing Einstein's general theory of relativity to de- 
scribe the evolution of the Universe, we base our uni- 
verse model on a relation between geometry and matter 
sources. A maximal reduction of this theoretical fun- 
dament is to consider the simplest conceivable geometry. 
Without putting in doubt that it might be an oversimpli- 
fication to assume a locally isotropic (and hence homoge- 
neous) geometry, standard cosmology conjectures the ex- 
istence of sources that would generate this simple geom- 
etry. As already remarked, the majority of these sources 
have yet unknown physical origin. Obviously, particle 
physicists take the demand for missing fundamental fields 
literally. But, as was emphasized above, the standard 
model has physical sense only, if a homogeneous-isotropic 
solution of Einstein's equations also describes the inho- 
mogeneous Universe effectively, i.e. on average. This is 
not obvious. The very fact that the distributions of mat- 
ter and geometry are inhomogeneous gives rise to back- 
reaction terms; we shall restrict them to those additional 
terms that influence the kinematics of the homogeneous- 
isotropic solutions. These terms can be viewed to arise 
on the geometrical side of Einstein's equations, but they 
may as well be put on the side of the sources. 

We start with a basic kinematical observation that lies 
at the heart of the backreaction problem. 



2. Non-commutativity 

Let us define spatial averaging of a scalar field \& on 
a compact spatial domain T> with volume Vt> '■= \T>\ 
through its Riemannian volume average 

{*{X\t)) v :=±- j ^(X\t)Ja?X ; J := y/det( 9ij ) . 

(6) 

The key property of inhomogeneity of the field \& is re- 
vealed by the commutation rule (24j |. (l4| : 



(7) 



where O denotes the trace of the expansion tensor; its av- 
erage describes the rate of volume change of a collection 



5 



of fluid elements along d/dt, 

d t Vv 



(©), 



Vd 



3H V 



(8) 



where we have introduced a volume Hubble rate H-p 
that reduces to Hubblc's function in the homogeneous 
case. Commutativity reflects the conjecture of the stan- 
dard model which implies that a realistically evolved in- 
homogeneous field will feature the same average charac- 
teristics as those predicted by the evolution of the (ho- 
mogeneous) average quantity; in other words, the right- 
hand-side of 10 is assumed to vanish. This rule also 
shows that backreaction terms deal with the sources of 
non- commutativity that are in general non-zero for in- 
homogeneous fields. 



3. Regional volume acceleration despite local deceleration 

Based on a first application of the above rule, we shall 
emphasize that there is not necessarily anti-gravity at 
work, e.g. in the 'redcapped' version of a positive cosmo- 
logical constant, in order to have sources that counteract 
gravity. Raychaudhuri's equation, if physically essential 
terms like vorticity, velocity dispersion, or pressure are 
retained, provides terms needed to oppose gravity, e.g., to 
support spiral galaxies (vorticity), elliptical galaxies (ve- 
locity dispersion) , and other stabilization mechanisms in- 
volving pressure (think of the hierarchy of stable states of 
stars until they collapse into a Black Hole). Admittedly, 
those terms are effectively 'small-scale-players'. Now, let 
us consider Raychaudhuri's equation restricted to irrota- 
tional dust, 



d t Q = A - 4ttG£> + 211 - I 2 



(9) 



with the principal scalar invariants of O a fc, 211 := 2/30 2 — 
2cr 2 and 1 := 6. Then, unless there is a positive cosmo- 
logical constant, there is no term that could counter- 
balance gravitational attraction and, at every point, 
dt@ < 0. Applying the commutation rule ([7|) for ^ = 0, 
we find that the averaged variables obey the same equa- 
tion as above despite non-commutativity [86] : 



d t (Q) v = A - AttG (g) v + 2 (11)^ - (I) 2 , 



(10) 



This result can be understood on the grounds that shrink- 
ing the domain I? to a point should produce the cor- 
responding local equation. Now, notwithstanding, the 
above equation contains a positive term that acts against 
gravity. This can be easily seen by rewriting the av- 
eraged principal invariants: we obtain 2 (II) .p — (I) 
2/3 ((6 - (6),,)% - 1/3 - 2 (a 2 ) v which 
pared with the corresponding local term 211 — I 2 = 
— 1/39 2 — 2a 2 , gave rise to a positive-definite averaged 
variance. It may appear 'magic' that the time-derivative 
of a (on some spatial domain T>) averaged expansion may 
be positive despite the time-derivative of the expansion 



com- 



at all points in T> is negative. As the above explicit calcu- 
lation shows, this is not 'unphysical' [5l[ but simply due 
to the fact that an average correlates the local contribu- 
tions, and it is this correlation (or fluctuation) that acts 
in the sense of a 'kincmatical pressure'. The interesting 
point is that this positive term is a 'large-scale player', as 
we shall make more precise below (the physical and obser- 
vational consequences of this term have been thoroughly 
explained and illustrated in the review paper pT6l| ) . 

What we should learn from this simple exercise is that 
any local argument, e.g. on the smallness of some per- 
turbation amplitude at a given point, is not enough to 
exclude regional ('global') physical effects that arise from 
averaging inhomogeneities. As we shall also learn below, 
such correlation effects must not be subdominant com- 
pared to the magnitude of the local fields, since they are 
related to the spatial variation of the local fields and, 
having said 'spatial', it could (and it will) imply a cou- 
pling to the geometry as a dynamical variable in Einstein 
gravitation. This latter remark will turn out very useful 
in understanding the potential relevance of backreaction 
effects in relativistic cosmology. 



4- The production of information in the Universe 

The above considerations on effective expansion 
properties can be essentially traced back to 'non- 
commutativity' of averaging and time-evolution, lying at 
the root of backreaction. (Note that additional backre- 
action terms that have been discussed in Subsection II B 21 
are also the result of a 'non-commutativity', this time 
between averaging and spatial rescaling - see also [1^].) 
The same reasoning applies to the 'production of infor- 
mation' in the following sense. Applying the commuta- 
tion rule ([7]) to the density field, W = g, 



(d t g) v - d t (g) v = 



d t S{g\\ (g) v } 
Vd 



(11) 



we derive, as a source of non-commutativity, the 
(for positive-definite density) positive-definite Lyapunov 
functional (known as Kullback-Leibler functional in infor- 
mation theory; [49( and references therein): 



s{e\\ (q) v } ■ = / gin 



Jd s X 



(12) 



This measure vanishes for Friedmannian cosmologies 
('zero structure'). It attains some positive time- 
dependent value otherwise. The source in (|1 lj) shows 
that relative entropy production and volume evolution 
are competing: commutativity can be reached, if the vol- 
ume expansion is faster than the production of informa- 
tion contained within the same volume. 

In [49( the following conjecture was advanced: 

The relative information entropy of a dust matter 
model S{g\ \ (g)^} is, for sufficiently large times, globally 
(i.e. averaged over the whole manifold £ that is assumed 
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simply-connected and without boundary) an increasing 
function of time. 



II. REFURNISHING THE COSMOLOGICAL 
EQUATIONS 



This conjecture already holds for linearized scalar per- 
turbations on a Friedmannian background (the growing- 
mode solution of the linear theory of gravitational insta- 
bility implies dt S > and S is, in general, time-convex, 
i.e. Of S > 0). Generally, information entropy is pro- 
duced, i.e. dt S > with 



d t S{g\\ (g) v } 
V v 



(6q6Q). 



(13) 

(and with the deviations of the local fields from their av- 
erage values, e.g. Sg := g — (g) v ), if the domain T> con- 
tains more expanding underdense and contracting over- 
dense regions than the opposite states contracting un- 
derdense and expanding overdense regions. The former 
states are clearly favoured in the course of evolution, as 
can be seen in simulations of large-scale structure. 

There are essentially three lessons relevant to the ori- 
gin of backreaction that can be learned here. First, struc- 
ture formation (or 'information' contained in structures) 
installs a positive-definite functional as a potential to in- 
crease the deviations from commutativity; it can there- 
fore not be statistically 'averaged away' (the same remark 
applies to the averaged variance of the expansion rate 
discussed before). Second, gravitational instability acts 
in the form of a negative feedback that enhances struc- 
ture (or 'information'), i.e. it favours contracting clusters 
and expanding voids. This tendency is opposite to the 
thermodynamical interpretation within a closed system 
where such a relative entropy would decrease and the sys- 
tem would tend to thermodynamical equilibrium. This 
is a result of the long-ranged nature of gravity: the sys- 
tem contained within T> must be treated as an open sys- 
tem. Third, backreaction is a genuinly non-equilibrium 
phenomenon, thus, opening this subject also to the lan- 
guage of non-equilibrium thermodynamics [72| , [7^ | , [85| . 
'Near-equilibrium' can only be maintained (not estab- 
lished) by a simultaneous strong volume expansion of the 
system. Later we discuss an example of a cosmos that 
is 'out-of-equilibrium', i.e. settled in a state far from 
a Friedmannian model that can be associated with the 
relative equilibrium state S = 0. 

In particular, we conclude that the standard model 
may be a good description for the averaged variables 
only when information entropy production is over- 
compensated by volume expansion, a property that is 
realized by linear perturbations on a FRW background. 
Thus, the question is whether this remains true in the 
nonlinear regime, where information production is ex- 
pected to be more efficient. 

Before we can go deeper into the problem of whether 
such backreaction terms, being well-motivated, are in- 
deed relevant in a quantitative sense, we have to study 
the governing equations. 



In this section we recall a set of averaged Einstein equa- 
tions together with alternative forms of these equations 
which put us in the position to study backreaction terms 
as additional sources to the standard Friedmann equa- 
tions. 



A. Averaged cosmological equations 

In order to find evolution equations for effective (i.e. 
spatially averaged) cosmological variables, we may put 
the following simple idea into practice. We observe 
that Friedmann's differential equations capture the scalar 
parts of Einstein's equations, while restricting them by 
the strong symmetry assumption of local isotropy. The 
resulting equations, Friedmann's expansion law (the en- 
ergy or Hamiltonian constraint) and Friedmann's accel- 
eration law (Raychaudhuri's equation), 



8irGg H - A 



3- + AnGgH - A 
a 



3/fc 
7a 



, 



(14) 



(15) 



can be replaced by their spatially averaged, general 
counterparts (for the details the reader is referred to 
[LH [TBI [ID, |2g|; an overdot abbreviates, in both cases, 
the covariant time-derivative d t = u^d^): 



a-D 
a-p 



lirG(g) 7 



A 



3— + AttG (g)^ - A 



Qv 



(16) 



(17) 



We have replaced the Friedmannian scale factor by the 
volume scale factor a-p, depending on content, shape and 
position of the domain of averaging 2?, defined via the 
domain's volume Vx>(t) — \T>\, and the initial volume 



1/3 



(18) 



Using a scale factor instead of the volume should not 
be confused with 'isotropy'. The above equations are 
general for the evolution of a mass-preserving, compact 
domain containing an irrotational continuum of dust, i.e. 
they provide a background-free and non-perturbative de- 
scription of inhomogeneous and anisotropic fields. (One 
could, of course, introduce an isotropic or anisotropic ref- 
erence background [24| or, explicitly isolate an averaged 
shear from the above equations to study deviations from 
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the kinematics of Bianchi-type models, as was done with 
some interesting results in [8|.) The new term appearing 
in these equations, the kinematical backreaction, arises as 
a result of expansion and shear fluctuations: 



and obtain: 



Q V := 2 (ll) v | (1)1 = 1 



(9 - (9), 



(19) 

I and II denote the principal scalar invariants of the ex- 
trinsic curvature, and the second equality follows by in- 
troducing the decomposition of the extrinsic curvature 
into the kinematical variables, as before. Also, it is 
not a surprise that the general averaged 3-Ricci curva- 
ture (1Z) -0 replaces the constant-curvature term in Fried- 
mann's equations. Note also that the term Qx> encoding 
the fluctuations has the particular structure of vanishing 
at a Friedmannian background, a property that it shares 
with gauge-invariant variables. 

In the Friedmannian case, Eq. (fl5|) arises as the time- 
derivative of Eq. (fT4|) . if the integrability condition of 
restmass-conservation is respected, i.e. the homoge- 
neous density qh oc a~ 3 . In the general case, however, 
restmass-conservation is not sufficient. In addition to the 
(built-in) integral 



(Q)v = 



Mr 



M v = M Vi , (20) 



we also have to respect the following curvature- 
fluctuation-coupling: 

-^d t (Q v a%) + ±-d t ({K) v al) =0. (21) 
a v a v 

This relation will be key to the understanding of how 
backreaction can take the role of Dark Energy. 



B. Alternative forms of the averaged equations 

We here provide three compact forms of the aver- 
aged equations introduced above, as well as some derived 
quantities. They will prove useful for our further discus- 
sion of the backreaction problem. 



8ttG (g) v - A = 



1 

4 ji> 



dt' Q 7 



dt 



al(t' 



(23) 

This equation is formally equivalent to its Newtonian 
counterpart [24[. It shows that, by eliminating the aver- 
aged scalar curvature, the whole history of the averaged 
kinematical fluctuations acts as a source of a generalized 
expansion law that features the 'Friedmannian part' on 
the left-hand-side of (|25|l . 



2. Effective Friedmannian framework 

We may also recast the general equations (|T6l fTTI) by 
appealing to the Friedmannian framework. This amounts 
to reinterpret geometrical terms, that arise through aver- 
aging, as effective sources within a Friedmannian setting. 
In the present case the averaged equations can be written 
as standard zero-curvature Friedmann equations for an 
effective perfect fluid energy momentum tensor with new 
effective sources 11511 : 



Q c « = (q)v ~ 



1 



l 



(K) 



a-D 



16ttG 16ttG 



16ttG" 



48ttG 



(24) 



8nGgf s - A 



3 2E+47rGG& + 3p£)-A 
a-D 



■ v 

f?eff 



3 {QeS 









(25) 



Eqs. |15]) correspond to the equations ([15]). fT7|) and (|2T|) . 
respectively. (Note that we may retain the constant- 
curvature part in the (scale-dependent) Friedmannian 
expansion law and instead consider only the deviations 
of averaged curvature from constant curvature as sources 
(see H|).) 



1. Generalized expansion law 

The correspondence between Friedmann's expansion 
law (TT4"]) and the general expansion law (flB)) can be made 
more explicit through formal integration of the integra- 
bility condition (|2Tj) : 

5^1 _ i_ f dt ' Qv ^(f) = 1 {{ K) V + Q v ) . 

a V a T> Jto al L 

The (domain-dependent) integration constant relates 
the new terms to the 'constant-curvature part'. We in- 
sert this latter integral back into the expansion law (TT^|) 



3. 'Morphed' Friedmann cosmologies 

In the above-introduced framework we distinguish the 
averaged matter source on the one hand, and averaged 
sources due to geometrical inhomogeneities stemming 
from extrinsic and intrinsic curvature (kinematical back- 
reaction terms) on the other. As shown above, the aver- 
aged equations can be written as standard Friedmann 
equations that are sourced by both. Thus, we have 
the choice to consider the averaged model as a (scale- 
dependent) 'standard model' with matter source evolving 
in a mean field of backreaction terms. This form of the 
equations is closest to the standard model of cosmology. 
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It is a 'morphed' Friedmann cosmology, sourced by mat- 
ter and 'morphed' by a (minimally coupled) scalar field, 
the morphon field [26[ . We write (recall that we have no 
matter pressure source here): 



with 



Q<s> 



2 v 



{0)v 



+ U V 



V 



V 



V 



(26) 



(27) 



where e = +1 for a standard scalar field (with positive 
kinetic energy), and e = — 1 for a phantom scalar field 
(with negative kinetic energy); (we have chosen the letter 
U for the potential to avoid confusion with the volume 
functional; e arises if we insist on a real-valued scalar 
field). Thus, in view of Eq. we obtain the following 
correspondence: 



1 



8ttG" 



- U v 



(28) 

Inserting ([28]) into the integrability condition (f2Tj) then 
implies that $£>, for $£> ^ 0, obeys the (scale-dependent) 
Klein-Gordon equation: 

d 

<f>v + 3H v <I> v + e——U{<f>v) = . (29) 

This correspondence allows us to interpret the kinemati- 
cal backreaction effects in terms of properties of scalar 
field cosmologies, notably quintessence or phantom- 
quintessence scenarii that are here routed back to models 
of inhomogeneities. For a full-scale discussion of this cor- 
respondence see [26j | . 



4- A note on closure assumptions 

This system of the averaged equations in the various 
forms introduced above does not close unless we specify a 
model for the inhomogeneities. Note that, if the system 
would close, this would mean that we solved the scalar 
parts of the GR equations in general by reducing them to 
a set of ordinary differential equations on arbitrary scales. 
Closure assumptions have been studied by prescribing a 
cosmic equation of state of the form p^ s = /3(g^ s ,aT>) 
[UGH, or by prescribing the backreaction terms through 
scaling solutions, e.g. Qt> oc a£> [26| . We shall come back 
to the important question of how to close the averaged 
equations in Subsect. IIII D 21 



C. Derived dimensionless quantities 

For any quantitative discussion it is important to pro- 
vide a set of dimensionless characteristics that arise from 
the above framework. 



1. The cosmic quartet 

We start by dividing the volume-averaged Hamilto- 
nian constraint (|16[) by the squared volume Hubble func- 
tional Hx> '■= a-p/a-D introduced before. Then, expressed 
through the following set of 'parameters' [87| . 

A 



8nG . , 

—F2\Q)V 



ZH. 



O v — - 



v 

ml 



3-fff, 



Q 6H 2 



(30) 



v 



the averaged Hamiltonian constraint assumes the form of 
a cosmic quartet [l6l [2l| : 



n 



V 



n 



= i 



(31) 



showing that the solution space of an averaged inhomo- 
geneous cosmology is three-dimensional in the present 
framework. In this set, the averaged scalar curvature pa- 
rameter and the kinematical backreaction parameter are 
directly expressed through (1Z) V and Qx>, respectively. 
In order to compare this pair of parameters with the 
'constant-curvature parameter' that is the only curva- 
ture contribution in standard cosmology to interpret ob- 
servational data, we can alternatively introduce the pair 



3of, 



(32) 



being related to the previous parameters by 



s 'fe ' lL QN 



9. 



n 



n 



V 



(33) 

After a little thought we see that both sides of this equal- 
ity would mimick a cosmological 'constant' in a Fried- 
mannian model. Note, in particular, that it is not the 
additional backreaction parameter alone that can play 
this role, but it is the joint action with the (total) cur- 
vature parameter, or, looking to the left-hand-side, it 
is the cumulative effect acquired during the history of 
the backreaction parameter. A positive cosmological con- 
stant would require this sum, or the effective history, re- 
spectively, to be positive. 



2. Volume state finders 

Like the volume scale factor a-p and the volume Hubble 
functional Hx>, we may introduce 'parameters' for higher 
derivatives of the volume scale factor, e.g. the volume 
deceleration functional 



T 2 ~ 2 m 



(34) 



aj) tL v 

Following [J, [73] (see also [HI and references therein) we 
may also define the following volume state finders involv- 
ing the third derivative of the volume scale factor: 



■d 



a-D 1 
a v Hi, 



n£(i+2ng)+2n2(i+4ng)-— o| 



(35) 
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and 



3(g c - 1/2) 



(36) 



The above definitions are identical to those given in 
[2> Iz3> however, note the following obvious and subtle 
differences. One of the obvious differences was already 
mentioned: while the usual state finders of a global ho- 
mogeneous state in the standard model of cosmology are 
the same for every scale, the volume state finders de- 
fined above are different for different scales. The other 
is the fact that the volume state finders apply to an in- 
homogeneous cosmology with arbitrary 3-metric, while 
the usual state finders are restricted to a FRW metric. 
Besides these there is a more subtle difference, namely a 
de gen eracy in the Dark Energy density parameter: while 
denote (with obvious adaptation) 1 — fi^ = Vl v 



degene 

SB 



v 



we have from Hamilton's constraint (|3 1 1) il x — aigT" K , 
i.e. so-called X-matter is composed of two physically 
distinct components. 



Cosmic equation of state and Dark Energy equation of 
state 



obeyed by the scale-dependent Friedmannian model. 
Again, a non-vanishing backreaction is associated with 
violation of 'equilibrium'. Note that a morphon field 
does not violate energy conditions as in the case of a 
fundamental scalar field (see [26|). Again it is worth em- 
phasizing that the above-defined equations of state are 
scale-dependent . 

With the help of these dimensionlcss parameters an in- 
homogeneous, anisotropic and scale-dependent state can 
be effectively characterized. 



III. IMPLICATIONS, FURTHER INSIGHTS 
AND PROSPECTS 



Having laid down a framework to characterize inho- 
mogeneous cosmologies and having understood the phys- 
ical nature of backreaction effects, does not entitle us to 
conclude about the quantitative importance of inhomo- 
geneities for the global properties of world models. It 
may well be that the robustness of the standard model 
also withstands this challenge. A good example is pro- 
vided by Newtonian cosmology that is our starting point 
for discussing the implications of the present framework. 



We already mentioned the possibility to characterize a 
solution of the averaged equations by a cosmic equation 
of state pf s = f3(gf s ,a v ) with wf s := pf s /gf s . Now, we 
may separate the morphon equation of state that plays 
the role of the Dark Energy equation of state 26] , 



x> = Qv - l/3(K)p 
W * ' Qv + (K) v 



(37) 



We can express the volume state finders through this 
equation of state parameter and its first time-derivative: 



m) 2H V ( 



and 



= 1 



v 



(39) 



As emphasized in 0, [77j], the above expressions have 
the advantage that one can immediately infer the case 
of a constant Dark Energy equation of state, so-called 
quiessence models, that here correspond to scaling solu- 
tions of the morphon field with a constant fraction of 
kinetic to potential energies [261 ] : 



-UvV v 



(11} i 



v 



(40) 



where the case Qv 
w% = -1/3 (i.e 
condition' 



(no kinematical backreaction) , or 
0$ +3pf = 0) corresponds to the 'virial 



2E kin 



pV 



(41) 



Thoughts on Newtonian cosmology and N— body 
simulations 



Analytical as well as numerical models for inhomo- 
geneities are commonly studied within Newtonian cos- 
mology. Essential cornerstones of our understanding of 
inhomogeneities rest on the Euclidean notion of space 
and corresponding Euclidean spatial averages. 



1. Global properties of Newtonian models 

The present framework can also be applied to the New- 
tonian equations and, indeed, at the beginning of its de- 
velopment the main result on global properties of New- 
tonian models was the confirmation of the FRW cosmol- 
ogy as a correct model describing the averaged inhomo- 
geneous variables. Technically, this result is due to the 
fact that the averaged principal invariants, encoded in 
Qv, are complete divergences on Euclidean space sec- 
tions and, therefore, have to vanish on some scale where 
we impose periodic boundary conditions on the deviation 
fields from the FRW background. The latter is a neces- 
sary requirement to obtain unique solutions for Newto- 
nian models (for details see [24T]). 

This point is interesting in itself, because researchers 
who have setup cosmological N-body simulations did not 
investigate backreaction: the vanishing of the averaged 
deviations from a FRW background was enforced by con- 
struction. The same remark applies to analytical mod- 
els, where a homogeneous background is introduced with 
the manifest implication of coinciding with the averaged 
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model, but without an explicit proof. The outcome that a 
FRW cosmology indeed describes the average of a general 
Newtonian cosmology correctly can be traced back to the 
(non-trivial) property that the second principal invari- 
ant II can indeed be written (like the first) as a complete 
divergence. Since this latter fact is not valid in Rieman- 
nian geometry, backreaction studies pose a clear task for 
generalizing current cosmological simulations and analyt- 
ical models. If backreaction is substantial, then current 
models must be considered as toy-models that have im- 
proved our understanding of structure formation, but are 
inapplicable in circumstances where the dynamics of ge- 
ometry is a relevant issue. We shall learn that these cir- 
cumstances are those needed to route Dark Energy back 
to inhomogeneities. 



B. Qualitative picture for backreaction— driven 
cosmologies 

If we just look at the definition of averaging, any value 
of a volume-averaged quantity is damped by volume ex- 
pansion. As for the backreaction term Q-p, the relevant 
positive term that could potentially drive an accelerated 
expansion in accord with recent indications from super- 
novae data [H [sll , is the averaged variance of the rate of 
expansion, cf. Eq. (|19p . This term, however, is quadratic 
and damped by the square of the inverse volume. How 
can we then expect that, in an expanding Universe, such 
a term can be of any relevance at the present time? Be- 
fore we give an answer to this question, let us make pre- 
cise the notion of 'relevance'. This can be done with the 
help of the averaged equations as has been first advocated 
by Kolb et al. (5a,H|. 



2. Backreaction views originating from Newtonian 
cosmology and relativistic perturbation theory of a FRW 
background 



As a rule of thumb we can say that any model (includ- 
ing relativistic models) that involves periodic boundary 
conditions imposed on fields on Euclidean space must 
be rejected as a potential candidate for a backreaction- 
driven cosmology. Such models are simply too restrictive 
to account for a non-vanishing (Hubble-scale) Q-p that 
is a generic property of relativistic models. Of course, 
also in those models, backreaction can be investigated (a 
detailed inve stig ation within Newtonian cosmology may 
be found in [25| as well as an application on the abun- 
dance statistics of collapsed objects |53[), but it is then 
only a rephrasing of the known cosmic variance within 
the standard model of cosmology. 

We refer to the term 'quasi-Newtonian' when we think 
of models that are restricted to sit close to a Friedman- 
nian state. Although we do not refer to the discussion of 
structure on super-Hubble scales [73j], 54 1, [6l|, [6{J, the 
following consideration would also apply there. The inte- 
grability condition (|2"Tj) , in essence, spells out the generic 
coupling of kinematical fluctuations to the evolution of 
the averaged scalar curvature. If that coupling is absent, 
this equation shows that Qt> oc VpJ 2 and (TZ) V oc a^ 2 , i.e. 
the averaged curvature evolves like a constant-curvature 
model, and backreaction decays more rapidly than the 



Acceleration condition 



averaged density, (g)^ oc . In other words, back- 
reaction cannot be relevant today (we shall make this 
more precise in the following). Therefore, as another rule 
of thumb, we may say that any (relativistic) model that 
evolves curvature at or in the vicinity of the constant- 
curvature model is rejected as a potential candidate for a 
backreaction-driven cosmology [26]. 

In summary, there is no Dark Energy that can be routed 
back to inhomogeneities on large scales in Newtonian and 
quasi-Newtonian models. 



Let us look at the general acceleration law ifTTj) . and 
ask when we would find volume acceleration on a given 
patch of the spatial hypersurface [55[ , [13, EH : 



3— = A - 4nG{g)v 



> 0. 



(42) 



We find that, if there is no cosmological constant, the 
necessary condition Q-p > AttG{q)t) must be satisfied on 
a sufficiently large scale, at least at the present time. 
This requires that Q-p is positive, i.e. shear fluctuations 
are superseded by expansion fluctuations and, what is 
crucial, that Q-p decays less rapidly than the averaged 
density [l7j]. It is not obvious that this latter condition 
could be met in view of our remarks above. We con- 
clude that backreaction has only a chance to be relevant 
(e.g. as defined through the inequality Eq. (|4"2")l today), 
if its decay rate substantially deviates from its 'Newto- 
nian' behavior and, more precisely, its decay rate must 
be weaker than that of the averaged density (or at least 
comparable, depending on initial data for the magnitude 
of Early Dark Energy [27[ , (Hi). 



2. Curvature-fluctuation coupling 

It is clear by now that a backreaction-driven cosmology 
[76| must make efficient use of the genuinly relativistic 
effect that couples kinematical fluctuations to the aver- 
aged scalar curvature, as is furnished by the integrability 
condition (f2"Tj) (or the Klein-Gordon equation (f2"9")l in the 
mean field description) . While the Newtonian and quasi- 
Newtonian frameworks suppress the scalar field degrees 
of freedom attributed to backreaction (or the 'morphon 
field' in the mean field description) by construction [5l| . 
and so cannot lead to an explanation of Dark Energy on 
the Hubble scale, general relativity offers a wider range 
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of possible cosmologies, since it is not constrained by the 
assumption of Euclidean or constant-curvature geometry 
and small deviations thereof. But, how can a cosmolog- 
ical model be driven away from a 'near-Friedmannian' 
state, if we do not already start with initial data away 
from a perturbed Friedmannian model? What is the 
mechanism of a backreaction-driven cosmology? 



3. The 'Newtonian anchor' 

Let us guide our thoughts by the following intu- 
itive picture. Integral properties of Newtonian and 
quasi-Newtonian models remain unchanged irrespective 
of whether fluctuations are absent or 'turned on'. Imag- 
ine a ship in a silent water and wind environment (ho- 
mogeneous equilibrium state). Newtonian and quasi- 
Newtonian models do not allow, by construction, that 
the ship would move away as soon as water and wind 
become more violent. This 'Newtonian achor' is lifted 
into the ship as soon as we allow for the coupling of fluc- 
tuations to the geometry of spatial hypersurfaces in the 
form of the averaged scalar curvature. It is this coupling 
that can potentially drive the ship away, i.e. change the 
integral properties of the cosmology. Before we are going 
to exemplify this coupling mechanism, e.g. by discussing 
exact solutions, let us add some understanding to the role 
played by the averaged scalar curvature. 

4- The role of curvature 

Looking at the integral of the curvature-fluctuation- 
coupling, Eq. (f!?2|) , we understand that the constant- 
curvature of the standard model is specified by the in- 
tegration constant k-r> K - This term does not play a cru- 
cial dynamical role as soon as backreaction is at work. 
Envisaging a cosmology that is driven by backreaction, 
we may as well dismiss this constant altogether. In such 
a case, the averaged curvature is dynamically ruled by 
the backreaction term and its history. Given this remark 
we must expect that the averaged scalar curvature may 
experience changes in the course of evolution (in terms 
of deviations from constant -curvature) , as soon as the 
structure formation process injects backreaction. (This 
picture is actually what one needs to solve the coinci- 
dence problem.) 

This mechanism can be qualitatively understood by 
studying scaling solutions, which impose a direct cou- 
pling, Qx> oc (7?.) p- (These scaling solutions correspond 
to quiessence fields, Eq. (|40p . and have been thoroughly 
studied by many people working on quintessence (see 
(78|, [35| and references therein.) In the language of a 
'morphon field', the mechanism perturbs the 'virial equi- 
librium', Eq. (|4Tj) . such that the potential energy stored 
in the averaged curvature is released and injected into 
an excess of kinetic energy (kinematical backreaction). 
In the relevant case of a 'Dark Energy scenario' the av- 



eraged curvature becomes less positive or more negative, 
respectively, for positive kinematical backreaction. Thus, 
in this picture, positive backreaction, capable of mimick- 
ing Dark Energy, is feeded by the global 'curvature energy 
reservoir'. It is clear that such a mechanism relies on an 
evolution of curvature that differs from the evolution of 
the constant-curvature part of the standard cosmology. 
Indeed, as we shall exemplify below, already a deviation 
term of the form {TZ) V — k-p.^^ 2 oc a^ 3 would be quanti- 
tatively sufficient to explain the needed amount of Dark 
Energy today. 

If we start with 'near-Friedmannian initial data', and 
no cosmological constant, then the averaged curvature 
must be negative today and of the order of the value that 
we would find for a void-dominated Universe [26| . Thus, 
curvature evolution is key to realizing backreaction ef- 
fects in the direction that we would need to explain Dark 
Energy. The difference to the concordance model is es- 
sentially that the averaged curvature changes from an 
almost negligible value at the CMB epoch to a cosmo- 
logically relevant negative curvature today. This is one 
of the direct hints to put backreaction onto the stage of 
observational cosmology. 

Let us add two remarks. First, it is not at all evi- 
dent that a flat Universe is necessarily favoured by the 
data throughout the evolution (Hoj . This latter analysis 
has been performed within the framework of the stan- 
dard model, and it is clear that in the wider framework 
discussed here, the problem of interpreting astronomi- 
cal data is more involved. Second, it is often said that 
spatial curvature can only be relevant near Black Holes 
and can therefore not be substantial. Here, one mis- 
takenly implies an astrophysical Black Hole, while the 
Schwarzschild radius corresponding to the matter content 
in a Hubble volume is of the order of the Hubble scale. In 
essence, a cosmologically relevant curvature contribution 
is tiny, but this property is shared by all cosmological 
sources. 

(The above qualitative picture is illustrated in detail in 
Rasanen's review f76j.) 

C. Exact solutions for kinematical backreaction 

The following families of exact solutions of the aver- 
aged equations are used to illustrate the mechanism of 
a backreaction-driven cosmology. Other implications of 
these examples are discussed in fl8j . 

1. A word on the cosmological principle 

We may separate the following classes of solutions into 
those solutions that respect the cosmological principle 
and those that do not. It is therefore worth recalling the 
assumptions behind the cosmological principle. There 
are various 'definitions' of this principle in the literature. 
The most realistic, however, refers to the existence of a 
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scale of homogeneity: in standard cosmology we assume 
that there exists a scale beyond which all observables do 
no longer depend on scale. It is beyond this scale where 
the standard model is supposed to describe the Universe 
on average; it is simply unreasonable to apply this model, 
even on average, to smaller scales, since the standard 
FRW model has an in-built scale-independence. On the 
same grounds, isotropy can only be expected on the ho- 
mogeneity scale and not below. Accepting the existence 
of this scale has strong implications, one of them being 
that cosmological parameters on that scale are represen- 
tative for the whole Universe. If this were not so, and 
generically we may think of, e.g. a decay of average char- 
acteristics with scale all the way to the diameter of the 
Universe, then also backreaction would differ outside the 
domain of homogeneity. That scale is thought to be well 
below the scale of the observable Universe and within our 
past-lightcone. Therefore, averaging over non-causally 
connected regions delivers the same values as those al- 
ready cumulated up to the homogeneity scale [22j. 

We are now briefly describing some exact solutions, 
and we mainly have in mind to learn about the coupling 
between curvature and fluctuations. 



2. Kinematical backreaction as a constant-curvature or a 
cosmological constant 

Kinematical backreaction terms can model a constant- 
curvature term as is already evident from the integrabil- 
ity condition (|21[) . Also, a cosmological constant need 
not be included into the cosmological equations, since 
Qt> can play this role [l6j], J2j|, [ja|, and even provides 
a constant exactly, as was shown in (55| and [la ]. The 
exact condition can be inferred from Eq. ((TTJ) and (f2"5|) 
and reads: 

4" I & Qv^aKt') = Q V , (43) 

a T> Jti at 

which implies Qx> = Qv(U) = const, as the only pos- 
sible solution. Such a 'cosmological constant' installs, 
however, via Eq. l|22j) . a non-vanishing averaged scalar 
curvature (even for k-p-, = 0): 

(TZ) V = ^ - 3Q V (U) . (44) 
a v 

This fact has interesting consequences for 'morphed' in- 
flationary models [58[. 



3. Out-of-equilibrium Einstein cosmos 

Following Einstein's thought to construct a globally 
static model, we may require the effective scale-factor as 
on a simply-connected 3-manifold £ without boundary 
to be constant on some time-interval, hence ds = ds = 



and Eqs. (JTTJ) and (fTBj) may be written in the form: 

Qs = 4ttG-^--A ; (K) s = 12ttG-^% + 3A , 

(45) 

with the global kinematical backreaction Qs, the globally 
averaged scalar 3-Ricci curvature (12) E , and the total 
restmass Ms contained in S. 

Let us now consider the case of a vanishing cosmolog- 
ical constant: A = 0. The averaged scalar curvature is, 
for a non-empty Universe, always positive, and the bal- 
ance conditions (|45p replace Einstein's balance conditions 
that determined the cosmological constant in the stan- 
dard homogeneous Einstein cosmos. A globally static 
inhomogeneous cosmos without a cosmological constant 
is conceivable and characterized by the cosmic equation 
of state: 

= 3Q S = const. => pf s = gf s = . (46) 

Eq. (|46p is a simple example of a strong coupling between 
curvature and fluctuations. Note that, in this cosmos, the 
effective Schwarzschild radius is larger than the radius of 
the Universe, 

a S = / S=r= F= = — 2GAfs = — aSchwarzschild j (47) 

hence confirms the cosmological relevance of curvature 
on the global scale £. The term 'out of equilibrium' 
refers to Subsection II C 4l in the above example volume 
expansion cannot compete with information production 
because the volume is static, while information is pro- 
duced (sec 18] for more details). 

Such examples of global restrictions imposed on the 
averaged equations do not refer to a specific inhomoge- 
neous metric, but should be thought of in the spirit of the 
virial theorem that also specifies integral properties but 
without a guarantee for the existence of inhomogeneous 
solutions that would satisfy this condition. 

4- A globally stationary inhomogeneous cosmos 

Suppose that the Universe indeed is hovering around 
a non-accelerating state on the largest scales. A wider 
class of models that balances the fluctuations and the av- 
eraged sources can be constructed by introducing globally 
stationary effective cosmologies: the vanishing of the sec- 
ond time-derivative of the scale-factor would only imply 
ds = const. =: C, i.e., as = &s + C(t — U), where the 
integration constant as is generically non-zero, e.g. the 
model may emerge [4(|, [4jJ from a globally static cos- 
mos, as '■= 1, or from a 'Big-Bang', if as is set to zero. 
In this respect this cosmos does not appear very different 
from the standard model, since it evolves at an effective 
Hubble rate £?s <xl/t. (There are, however, substantial 
differences in the evolution of cosmological parameters, 
see [3.) 



13 



The averaged equations deliver a dynamical coupling 
relation between and (72.) s as a special case of the 
integrability condition (|2Tj) : 



1 4C 3 
d t <2 s + -d t (7e) s - — 



(48) 



The cosmic equation of state of the A— free stationary 
cosmos and its solutions read [Fa, E3l : 



1 



(49) 
(50) 



The total kinematical backreaction QsVs = 47rGMs is 
a conserved quantity in this case. 

The stationary state tends to the static state only in 
the sense that, e.g. in the case of an expanding cosmos, 
the rate of expansion slows down, but the steady increase 
of the scale factor allows for a global change of the sign 
of the averaged scalar curvature. As Eq. (|50[) shows, 
an initially positive averaged scalar curvature would de- 
crease, and eventually would become negative as a result 
of backreaction. This may not necessarily be regarded as 
a signature of a global topology change, as a correspond- 
ing sign change in a Friedmannian model would suggest 
(see Subsect: lmTJT|) . 

The above two examples of globally non-accelerating 
universe models evidently violate the cosmological prin- 
ciple, while they would imply a straightforward explana- 
tion of Dark Energy on regional (Hubble) scales: in the 
latter example the averaged scalar curvature has acquired 
a piece oc a^ 3 that, astonishingly, had a large impact on 
the backreaction parameter, changing its decay rate from 



oc a 



s to oc o s . This is certainly enough to produce 



sufficient 'Dark Energy' on some regional patch due to 
the presence of strong regional fluctuations. However, 
solutions that respect the cosmological principle and, at 
the same time, satisfy observational constraints can be 
constructed. As an example, scaling solutions [2(| have 
been exploited for such a more conservative approach, 
also outlined in Subsect. IIII Dl 



5. Explicit inhomogeneous solutions 

If we wish to specify the evolution of averaged quanti- 
ties without resorting to phenomenological assumptions 
on the equations of state of the various ingredients, or on 
the necessarily qualitative analysis of scaling solutions, 
or with specific global assumptions, we have to spec- 
ify the inhomogeneous metric. Natural first candidates 
are the spherically-symmetric Lemaitre-Tolman-Bondi 
(LTB) solutions that were first employed in the context 
of backreaction in [H| and [74[ • 

Considerable effort has been spent on LTB solutions 
and, especially recently, relations to integral properties 



of averaged cosmologies have been sought. Interestingly, 
[H[ also found a strong coupling between averaged scalar 
curvature and kinematical backreaction, and LTB solu- 
tions also feature an additional curvature piece oc a^ 3 
on some domain T>. There are obvious shortcomings of 
LTB solution studies that consider the class of on aver- 
age vanishing scalar curvature, since in that class also 
Qt> = (67[; also here, a non-vanis hing averaged curva- 
ture is crucial to study backreaction [30] . However, there 
is enough motivation to study the extra effect of aposi- 
tive expansion variance to fit observational data ([33j and 
references below). 

The value of LTB studies is more to be seen in the 
specification of observational properties such as the lu- 
minosity distance in an inhomogeneous metric 0, H, H 



12], |43|, |45| . Although interesting results were obtained, 



especially in connection with the interpretation of su- 
pernova data, care must be taken when determining e.g. 
just luminosity distances, since the free LTB functions 
may fit any data (62J. Generally, apart from mistakes 
(e.g. setting the shear to zero), those studies sometimes 
confuse integral properties of a cosmological model with 
local properties (e.g. the scale factor aj> and a local scale 
factor in the given metric form). The averaged equations 
cannot predict luminosity distances unless one considers 
averages on the lightcone (see, however, the strategy by 
Palle [6a|), which in turn is related to the issue of light 
propagation in an inhomogeneous Universe (see [52| , [lll | , 



and discussion and references in 



])• 



D. Where to go from here? 

In this subsection we are going to outline several strate- 
gies towards the goal of understanding the quantitative 
importance of backreaction effects, and to device meth- 
ods of their observational interpretation. All the topics 
listed are the subject of work in progress. 



1. Global aspects 

The question of what actually determines the aver- 
aged scalar curvature is an open problem. For a two- 
dimensional Riemannian manifold this question is an- 
swered through the Gauss-Bonnet theorem: the av- 
eraged scalar curvature is determined by the Euler- 
characteristic of the manifold. Hence, it is a global topo- 
logical property rather than a certain restriction on local 
properties of fluctuations that determines the averaged 
scalar curvature. If such an argument would carry over 
to a three-dimensional manifold, then any local argu- 
ment for an estimate of backreaction would obviously be 
off the table. (There are related thoughts and results 
in string theory that could be very helpful here.) In 
ongoing work [22, Hl| we consider the consequences of 
Perelman's proofs that were mentioned in Subsect. II B 21 
There is no such theorem like that of Gauss and Bonnet 
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in 3D, but there are uniformization theorems that could 
provide similar conclusions. For example, for closed in- 
homogeneous universe models we can apply Poincare's 
conjecture (now proven by Perelman) that any simply- 
connected three dimensional Ricmannian manifold with- 
out boundary is a homcomorph of a 3-sphere. Ongoing 
work concentrates on the multi-scale analysis of the cur- 
vature distribution and the related distribution of kine- 
matical backreaction on cosmological hypersurfaces that 
feature the phenomenology we observe. All these stud- 
ies underline the relevance of topological issues for a full 
understanding of backreaction in relativistic cosmology. 
To keep track with the developments in Riemannian ge- 
ometry and related mathematical fields will be key to 
advance cosmological research. In this line it should be 
stressed that the averaged scalar curvature is only a weak 
descriptor for the topology in the general 3D case, and 
information on the sectional curvatures or the full Ricci 
tensor is required. In observational cosmology there are 
already a number of efforts, e.g. related to the observa- 
tion of the topological structure of the Universe derived 
from CMB maps (for further discussion see Jl8| and for 
topology-related issues see [H, [H| , [36| , [j| \(^). 



2. Approximation methods for backreaction 

There is a large body of possibilities to construct a 
generic inhomogeneous metric. First, there is the pos- 
sibility of using standard methods of perturbation the- 
ory. Although the equations and 'parameters' discussed 
in this lecture can live without introducing a back- 
ground spacetime, a concrete model for the backreac- 
tion terms can be obtained by employing perturbation 
theory (preferably of the Lagrangian type) and, hence, 
a reference background must be introduced. But, the 
construction idea is (i) to only model the fluctuations by 
perturbation theory ( the term Qr> ) and to find the final 
( non-perturbative ) model by employing the exact frame- 
work of the averaged equations. Such a model is currently 
investigated by paraphrasing the corresponding Newto- 
nian approximation [2 51 ] . Second, we could aim at finding 
an approximate evolution equation for Qx> by (ii) clos- 
ing the hierarchy of ordinary differential equations that 
involve the evolution of shear and the electric and mag- 
netic parts of the projected Weyl tensor. In this line, (hi) 
further studies of cosmic equations of state (like, e.g., the 
Chaplygin state [ill) are not only a clearcut way to close 
the averaged equations, but also a way to classify differ- 
ent solution sectors. All these models could be subjected 
to (iv) standard dynamical system 's analysis to show their 
stability in the phase space of their parameters [80L [83[ . 
Note in this context that the FRW cosmology as an av- 
eraged model is found to be stable in many cases, but 
there is an unstable sector that just lies in the right cor- 
ner needed to explain 'Dark Energy', as was found in 



3. Issues of interpretation of backreaction within 
observational cosmology 

Here, the most important step will be (i) to investigate 
the averaging formalism on the lightcone. Such a frame- 
work is currently constructed. It relates not only to all 
aspects of observations in terms of distances within inho- 
mogeneous cosmologies, but also links directly to initial 
data in the form of, e.g., CMB fluctuation amplitudes and 
the integrated Sachs- Wolfe effect. Relating lightcone av- 
erages to cosmological model averages is also possible and 
is in the focus of this investigation. The consequences of 
a quantitative importance of an integrated backreaction 
history on the lightcone are obvious. Applying generic 
redshift -distance relations e.g. to galaxy surveys would 
put us in the position to better understand the actual 
distribution of galaxies that are currently mapped with 
the help of FRW distances. If expansion fluctuations are 
a dominant player on large-scales, we can imagine that 
also the galaxy density maps would be affected. 

Several times we have already pointed out that (ii) 
scale-dependence of observables is key to understand the 
cosmological parameters in the present framework. View- 
ing observational data with this additional discrimination 
power of a scale-dependent interpretation of backreac- 
tion effects, there is furthermore (iii) a link to the Dark 
Matter problem that certainly is an important task to 
be formulated, and some more words on this particular 
problem are in order. 



4- A common origin of Dark Energy and Dark Matter? 

Concentrating on the Dark Energy problem has led us 
to focussing on a positive contribution of Qt> on large 
scales. However, the kinematical backreaction Qx> it- 
self can also be negative, and a sign-change may actu- 
ally happen by going to smaller scales. Looking at the 
phenomenology of large-scale structure reveals strongly 
anisotropic patterns, so that it is not implausible that 
on the scales of superclusters of galaxies we would find a 
shear-dominated Qr> < (this was actually found in the 
Newtonian investigation 25] that, however, suffers from 
the fact that Qx> is restricted to drop to zero on the peri- 
odicity scale of the fluctuations). Thus, again as a result 
of its scale-dependence, the kinematical backreaction pa- 
rameter can potentially be the origin of a Kinematical 
Dark Energy, but also of a Kinematical Dark Matter [l3T ] . 

Mapping kinematical backreaction with a 'morphon 
field' opens further links to previous studies that tried 
to model Dark Energy and Dark Matter by a scalar field 
( 65] and references to earlier work therein). 

With this in mind, the volume deceleration functional 
([M]) can change sign too, but this crucially depends on 
the value of the matter density parameter fi^. Let us 
now discriminate different spatial scales. We denote by 
7i the homogeneity scale, by £ the scale of a typical void, 
and by M the scale of a matter-dominated region (e.g. 
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the scale of galaxy clusters) [22| . We notice from Eq. (|34|) 
that, for a small value of fi^, a smaller negative value of 
fig is needed to obtain volume acceleration on the ho- 
mogeneity scale, q 1 ^ < 0. The precise criterion for vol- 
ume acceleration today is — f2g° > 0^°/4. Thus, for 
fl™ » 1/4 we would need -fig > 1/16, while for a 
model with no (exotic) Dark Matter, e.g. fij; k 1/30 
we would just need — fig > 1/120. This is a game with 
numbers until we understand the scale-dependence of 
the matter distribution. We know that on smaller scales 
(from galaxy halos to clusters and superclusters of galax- 
ies) we would need fi^° ~ ~ 1/4, and here there is 
certainly a problem to explain this amount by baryonic 
dark matter only. At the moment we can just mention 
this problem and say that a quantitative investigation of 
scale-dependence of the density parameter is an impor- 
tant task: if no non-baryonic Dark Matter is assumed, it 
would require to reconcile a rather large value on smaller 
scales with a small value on large scales. In a forthcom- 
ing work this issue is further discussed and investigated 
|22 |; it is indeed the case that the matter density pa- 
rameter drops substantially at around the scale £ in a 
cosmological slice that is volume-dominated by voids. 



picture, but ii) a quantitative problem to reconcile this 
picture with the data in the sense that there is not enough 
time for the mechanism to be sufficient. In that situation 
we 'lost' the standard model for a correct description of 
the late Universe, and we do not reach an explanation, 
unless we allow for initial data that are non-standard. 
This in turn asks for a comprehensive understanding of 
these required initial data, hence reconsideration of in- 
homogeneous inflationary models and their fluctuation 
spectrum at the exit epoch. We opened Pandora's Jar. 

Notwithstanding, I would consider such a situation as 
the beginning of a fruitful development of cosmology. 
As previously mentioned, the issues of scale-dependence 
of observables, the priors underlying interpretations of 
observations, the large subject of Dark Matter and, of 
course, the issue of Dark Energy, will be all interlocked 
and ask for a comprehensive realistic treatment beyond 
crude idealizations. 

Even if we would 'only' find a 10 percent effect, rather 
than 75 percent, these studies would have justified their 
quantitative importance for observational cosmology, and 
what is to be expected, would substantially improve our 
understanding of the Universe. 



E. A short conclusion: opening Pandora's Jar 

Let us conclude by stressing the most important issue: 
quantitative relevance of backreaction effects. Especially 
the recent efforts, spent on the backreaction problem by 
a fairly large number of researchers, added substantial 
qualitative understanding to the numerious previous ef- 
forts that were undertaken since George Ellis initiated 
this discussion in 1984 |38( (see references in 39]). The 
issue remains unresolved to date: an explanation of Dark 
Energy along these lines is attractive and also physically 
plausible, but a reliable and unambiguous estimate of the 
actual influence of these effects is lacking. This situation 
may change soon and for this to happen it requires con- 
siderable efforts, for which some possible strategies have 
been outlined in the last subsection. 

After those results are coming in, we may face a more 
challenging situation than anticipated by the qualitative 
understanding that we have. For example, while the 
explanation of Dark Energy by quintessence (or phan- 
tom quintessence) still allows to hide the physical conse- 
quences behind a scalar field that is open for a number of 
explanations, the mapping of a scalar field to the backre- 
action problem, as outlined in Subsection III B 3[ can no 
longer keep a phenomenological status: fluctuations ex- 
ist and can be measured. There are no free parameters, 
there are initial data that can be constrained. Despite 
being premature, it is nevertheless allowed to speculate 
that the outcome is i) a confirmation of the qualitative 
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